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Introduction 


The purpose of this report is to inform NASA Headquarters and 
NASA- Langley on the work completed on grant NSG-1046 entitled, "An 
Analytical Study of Effects of Aeroelasticity on Control Effectiveness", 
during the period June 1974 through March 19n. This is the final 
report submitted to NASA on this grant. 

Completed Work and Conclusions 

The description and conclusion of the completed work is divided 
into ten small studies. 

1 . Effect of Changing Constant Percent Chordwise Lines on ACp Values 

This study was conducted to find the best paneling scheme which 

could be used in the Elastic Ability Derivative (ELASTAD) program 

of Reference 1. This program predicts AC values over two and three 

dimensional elastic wings, with and without camber and flap, at subsonic 

and supersonic speeds. The program of Reference 1, based on the theory 

of Reference 2, uses equidistant or almost equidistant constant 

percent streamwise lines (CPSWL's). Therefore, this study considers 

only the effect of changing constant percent chordwise lines (CPCWL's) 

on AC„ values. 

P 

Several different chordwise paneling schemes v-jere used, on two 
and three dimensional wings with and without flaps, to predict ACp 
distributions using ELASTAD program. These ACp distributions were 
compared against analytical and/or experimental distributions. The 
best AC distributions were obtainea by using the Modified Woodward 
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Scheme for constant percent chordwise lines with the control point 
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located at "0.85" of the local panel chord. In the (<tod1f1ed Uoo<kard 
Scheme the panels nearest the leading and trailing edges are half the 
size of the panels in-between. For wings witK flaps* the Modified 
Woodward Scheme should be used separately* ahead and behind the flap 
hinge line. A Cp distributions over thick airfoils, like the 
6A(W)-1* should be used with caution because at non-zero angles of 
attack the predicted ACp values are higher at the trailing edge than 
those obtained by experiment. At zero angle of attack, the difference 
between predicted and experimental A Cp values is significant both at 
the leading and the trailing edge. 

The details of this study are given in reference 3. Further work 
on control point location is reported in section 9. 

2. Structural Complexity Study for Evaluation of Structural Influence 
Coefficien t Matrices 

Several wings of either solid or built-up type construction were 
used for this study, the details of which are given in reference 4. 

The structurc.1 influence coefficient matrices for these v-zings were 
evaluated by using ELASTAD- and NASTRAM-program (ref. 5). ELASTAD 
program uses slender beam theory to represent the wing, whereas NA5TRAN 
uses the actual elastic properties of the wing. 

In ELASTAD program the wing is represented by an elastic axis, 
thus the input data needed for a wing is* 

a. Unit loading point locations. 

b. Elastic axis coordinates. 

c. Bending stiffness (El) of elastic axis segments. 
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d. Torsional stiffness (6J) of elastic axis se^nts. 

e. lASIGN-array, which assigns aerodynamics panels (loading 
points) to elastic axis end points. 

A computer program was written to provide elastic axis coordinates and 
El- and GJ-values for its segments for solid wings (see Appendix in 
ref. 4). The unit loading point locations and lASIGN-array was 
evaluated by hand for solid wings. For built-up wings, all the input 
data was hand calculated. 

In NASTRAN program, the solid wings were assumed to be con^iosed 
of several triangular and/or quadrilateral plate element. An average 
thickness for each element was calculated by dividing its volume by 
its planform area. The built-up wings were modeled by using triangular 
and/or quadrilateral plates, shear panels and rods. The upper and 
lower surfaces were represented by plates, the spars and ribs by shear 
panels, and the flanges by rods. 

Most of the experimental results presented deflection influence 
coefficients (OIC's) only. The FLASTAD program, which uses slender 
beam theory, calculates streamwise rotational influence coefficients 
(RIC's); whereas, the NASTRAN program, which employs the actual elastic 
properties of the wing, calculates both OIC's and RIC's. The appropriate 
comparisons for all the wings are shown in reference 4. 



4 


In case of solid wings, the DIC's obtained by using NASTRAM program 
matched reasonably well with experimental measurements. In cases of 
built-up wings, the NASTRAN predicted slightly lower deflections than 
the experiment. These larger experimental deflections could be due 
to the mechanical construction of the wings. Rotational influence 
coefficients obtained by using ELASTAD and NASTRAN programs do not 
match exactly at any particular point on all the wings studied. The 
smallest difference between the two sets of rotations occurs, for 
untapered solid- and built-up - wings of aspect ratio between 2 and 
6, when load- and rotation-points are on the elastic axis. For the 
load Or rotation-points away from the elastic axis, this difference 
increases but is still reasonable. The larger difference is due to 
the SBM assuming the rigid links between the elastic axis endpoints 
and the load- and rotation-points; whereas, the NASTRAN employs the 
actual elastic properties of the structure. Thus, it is suggested to 
use the NASTRAN program for evaluating the RIC's to be used in 
ELASTAD program. 

3. Arrow V/inq Study 

H. W. Carlson concluded from a study of arrow wings (ref. 6) that 
the discrepancy betv^een the experimental and linearized theory estimates 
of pressure distribution might be due to aerolastic deflections and 
the presence of vortex flow. An arrow wing of symmetric airfoil 
section of reference 6 was studied for its aeroelastic effects (see 
ref. 7). 

The structural influence coefficient matrix for this v/ing was 
calculated by using the NASTRAN program (ref. 5). This was done by 
dividing the wing into triangular- and quadrilateral -plate elements. 
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The NteTBAN program can only handle constant thickness plate elements 
and so an average thickness of each elaient is used for calculating the 
structural influence coefficient matrix. This matrix was manipulated to 
generate another matrix which conformed to the aerodynamic paneling. 

Mass of each of the panel was calculated by multiplying its area to mean 
thickness and density. 

All the above mentioned information was used in ELASTAD program to 
calculate rigid and elastic A Cp distributions at the experimental angle 
of attack (8 degrees) and Mach number (2.05). The difference between the 
rigid and elastic A Cp's was small, which meant elastic effects were 
small. These A Cp's were quite different from those measured experi- 
mentally. 

It was of interest to compare the streanrwise rotations of panels 
obtained theoretically, by using ELASTAD program, and those obtained by 
employing experimental A Cp loading. The difference between the two 
sets of rotations is small at the root and large at the tip. The largest 
rotation of 1.94 degrees was predicted by experimental loading and 3.06 
degrees by ELASTAD program. Both experimental and theoretical rotations 
are small which means elastic effects are small. Thus, the discrepancy 
between the experimental ar.d theoretical Cp values can be attributed 
to the vortex flow. 

4. Modification of ELASTAD Program 

Geometry part (AEREAD) of ELASTAD program has been modified to plot 
top and side views and structural part (AERELAS) to plot panels and elastic 
axis of the models. Until recently lASIGN - array had lo be done graphically 
and punched for use in ELASTAD. A subroutine hos been written and tested 
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to do this automatically. This makes input data preparation easier and 
panel assignment consistent. 

5. Structural Matrix Conversion Package 

A computer program was written to convert a structural influence 
coefficient matrix from its structural network to the one which conforms 
to the aerodynamic paneling of ELASTAD programi. The details of this pro- 
gram are given in reference 8. This program was used in Transonic 
Aircraft Technology project, which is described next. 

6. Participation in the Transonic Aircraft Technology (TACT) Project 
In this project the rigid and elastic stability derivatives were 

caldulated for TACT aircraft for the following flight conditions: 

H = 0.60, = 300, 500, 533, 600, 825, 1050 Ibs/ft^ 

M = 0.85, = 300, 600, 825, 1050 Ibs/ft^ 

M = 0.90, q^ = 300, 600, 825, 1050 Ibs/ft^ 

All the stability derivatives are listed and plotted in reference 9. 

The following observations are noted from the results. 

The experimental rigid lift curve slope (C. ) is always higher than 

'a 

that calculated by ELASTAD and the difference betvs/een the two increases 

with Mach number. At M = 0.6, the difference between the experimental 

and ELASTAD values is 5 percent which is assumed to be in the error 

bound, but at M = 0.85, the difference increases to 10 percent. The 

difference between experimental and ELASTAD numbers could be due to 

unmodeled transonic effects being significant and possibly the vortex 

flow due to strake type behavior of the planform. At a constant Mach 

number elastic and lift coefficient for zero angle of attack ) 
a 0 
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decrtftse pressure* This is due to el as tic unloading of 

the aircraft for svmpt-back wings and results from the fact that the 
Increase In local a associated with twisting Is overcome by the decrease 
In local o due to bending. Ihe non-zero mass C. 1$ higher than zero 

mass C. . These higher magnitudes of C. Indicates that the flap region 

of the wing and tall have larger mass concentration and so the Increase 
In local ot due to twist Is higher than decrease In local a due to bending 
as conqsared to zero mass . This could also be explain^ by examining 


the variation of lift coefficient with load factor 


which are 


positive. A positive value of indicated that the lift coefficient 

an" 

Increases as the load factor increases, which, for this wing, could happen 
only when the flap region and tail are heavier. When a positive load 
factor is applied to such an aircraft, the local a increases due to the 
inertial forces which act opposite to the direction of motion. Thus, 
the lift is also increased. The equation for non-zero mass C. 

ct 

HP 

(Ref. 9) suggests that for a positive value of L , the non-zero mass 

dn 

C. is always higher than zero mass C. . 

S Ss 

The experimental rigid pitching nwinent curve slope ) is always 

a 

less negative than the one calculated by ELASTAD and the difference 
between the two increases with Mach number. At M = 0.6, the difference 
between the experimental and ELASTAD values is 3 percent, but at M ^ 0.85 
the difference is 60 percent. The difference between experimental and 
ELASTAD numbers could be due to the same unmodeled effects previously 
mentioned. At a constant Mach number, the zero-mass elastic becomes 
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less negative as the (lynamic pressure Increases and the pitching nonient 

coefficient for zero angle of attack ) is less negative than the 

"o 

rigid value at low dynamic pressure, but becomes n»re negative at higher 
dynamic pressures. These variations are again due to elastic unloading 

of the aircraft with swept-back wings. The non-zero mass is always 

a 

nore negative than zero mass . The same logic of flap region and 

tail having larger concentration of mass applies liere also. The variable 
to be noted here is the variation of pitching moment coefficient with 

load factor ^^m , which is negative. A negative value for *^^m 
dn dn 

always means that non-zero mass will be irore negative than zero 

a 

mass C_ , unless is larger than the lift coefficient for trimning 
""a dn 

HP 

which is highly unlikely because is always a small number. 

dn 

The experimental rigid static margin (C^ ) is always less 

a -a 

negative stable than the one calculated by ELASTAD and the difference 
between the two increases with Mach number. At M = 0.6, the difference 
between the experimental and ELASTAD values is 10 percent, but at M ■ 0.85 
the difference is 67 percent. The difference between experirental and 
ELASTAD numbers could be due to the same unmodeled effects previously 
mentioned. At a constant Mach number, the zero-mass elastic static 
margin becomes less negative as the dynamic pressure increases. The 
non-zero mass static margin is always more negative than zero mass static 
margin indicating that the effect of masses is to stabilize the aircraft. 
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7. Ex^nsioff of Stfuctoral Influenca Cotfflcltnt tutrix Evaluatton 
Packigt to Siwot-fonwinl W1n« 

Tmo trrors In tha ELASTM program mr^ detected during a routine 
calculation of a stroctural Influence coefficient natrlx. These »«ro; 

a. Wrong coding of few Fortran statements 

b. Wrong definition of one of the Input variables 

After making these nodlflcatlons. It was roallzed that ELASTAD sNiuld 
te modified further te Include the evaluation of structural Influoice 
coefficient iratrlx for swept forward wings and for the cases when 

f 

panel controlds on the horizontal tall lie behind the last panel centroid 
of the fuselage. All these modifications have been checked out. 

8. Leading-Edge Vortex Separation Study 

A nimdierlcal method Is developed to predict distributed and total 
aerodynamic characteristics for low aspect-ratio wings with partial 
leading-edge separation. The flow Is assumed to be steady and Inviscld. 
The wing boundary condition Is fornwlated by the Quasl-Vortex-Lattlce 
n»thod. The leading-edge separated vortices are represented by discrete 
free vortex elements which are aligned with the local velocity vector 
at mid-points to satisfy the force free condition. The wake behind the 
trail ing-edge Is also force free. The flow tangency boundary condition 
Is satisfied on the wing. Including the leading-and trail ing-edges. 
Comparison of the predicted results with complete leading-edge separation 
has shown reasonably good agreement. For cases with partial leading-edge 
separation the lift Is found to be highly nonlinear with angle of attack. 

The theoretical details of this study are given In reference 10 and 
the computer program In roference 11. This program was recently described 
In the Fall 1979 Issue of the quarterly publication "NASA Tech Briefs”. 
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9. Evaluation of HIng-TIp-Suctlon at Subsonic and Supersonic Speeds 

Tho aerodynamic method which has been used In EL/^TAD program Is 

generally known as Woodward's panel nmthod. A simplified version of 
ELAST/U) program was modified to check out the concept of evaluating the 
leading-edge and side-edge suction forces. 

Woodward's panel method for subsonic and supersonic flow Is In^roved 
by e»i^loy1ng control points determined by exactly matching two-dimensional 
pressure at a finite number of points. The results show great Improvement 
in the predicted pressure distribution of a flapped airfoil. With the 
paneling sche«« of cosine law in both chordwise and spanwise directions, 
the method Is shown to accurately predict leading-edge and side-edge 
suction forces of various configurations In subsonic and supersonic flow. 

Based on the extensive comparison of present prediction with other 
theoretical results, It may be concluded that the present Improved 
Woodward's panel n^thod is generally accurate in predicting the leading-edge 
and side-edge suction forces and the centers of these forces In subsonic 
and supersonic flow. The good accuracy of the present method has also 
been demonstrated for cairijered and flapped airfoils. Because of generality 
of the panel method, the present improved method can therefore be used 
not only to predict the vortex Htt of complex planforms through the 
method of suction analogy, but also to calculate certain lateral -directional 
stability derivatives as well. 

The details of this study are given in references 12 and 13. 

1 0 , Evaluation of Peflection Influc .-.e Coefficient Matrix for Solid 
Wing Model 

Several solid vn*ng nwdels, for which deflection data was available, 
were initially used for this study. Later, on a solid wing of transport 
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type aircraft was used to measure deflection for several load conditions. 
These wings were modeled on Structural Performance Analysis and Redesign 
(SPAR) program (ref. 14) by using first plate el«nents and then solid 
elements (ref. 15 and 16). The solid element representation of the wings 
resulted in a better correlation of measured and calculated deflections. 

By using solid el®i»nt representation the deflection Influence coefficient 
matrices can be calculated with about errar. 
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